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The energy conservation is a general law of nature. In the classical physics, the energy Wab 
of a conservative system AB that contains the objects A and B is equal to a sum of the positive 
energies Wa and Wb of the isolated objects A and B, Wab = Wa + Wb- We show that the 
energy conservation does not exhibit the "classic law" {Wab 7^ Wa + Wb) if the physical objects 
are waves or they do have a wave nature of microscopic (quantum) objects. The "extra energy", 
" negative energy" and " annihilation of energ'ij^ are predicted for multi-wave (multi-beam) systems. 
The paradoxical phenomenon is demonstrated in context of the extraordinary transmission of light 
and matter through subwavelength aperture arrays assisted by surface waves [T.W. Ebbesen et al., 
Nature (London) 391, 667 (1998)) and E. Moreno et al., Phys. Rev. Lett. 95, 170406 (2005)]. The 
physical mechanism of "non-classic energy conservation" provides the transmission enhancement 
(" extra energy" ) not by coupling to the surface waves but by interference between diffracted waves 
(beams) in the far field. Furthermore, the " annihilation of energy" by destructive interference of 
fields at the detector leads to the zero transmission. Another surprising prediction of the model is a 
weak Wood's anomaly, which is present in a classical Young type two-slit system. The Wood anoma- 
lies in transmission spectra of optical gratings, a long standing problem in optics, follows naturally 
from interference properties of the model. We believe that the "non-classic energy conservation" 
demonstrated in the present paper gains physical inside into the nature of energy conservation. The 
physical mechanism of the " extra energy" , " negative energy" and " annihilation of energy" affects the 
energy conservation in all subwavelength systems, for an example, in nanooptics and nanophotonics 
based on transmission and beaming of light or matter by subwavelengths apertures. 
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The scattering of waves by apertures is one of the ba- 
sic phenomena in the wave physics. The most remarkable 
feature of the light scattering by subwavelength apertures 
in a metal screen is enhancement of the light by excitation 
of electron waves in the metal. Since the observation of 
enhanced transmission of light through a 2D array of sub- 
wavelength metal nanoholes Q , the phenomenon attracts 
increasing interest of researchers because of its poten- 
tial for applications in nanooptics and nanophotonics . 
Recently, the enhanced transmission through subwave- 
length apertures was predicted also for matter waves Q . 
The enhancement of light is a process that can include 
resonant excitation and interference of surface plasmons 
H S H H , Fabry-Perot-like intraslit modes 0, i, H [13 > 
and evanescent electromagnetic waves at the metal sur- 
face ■ In the case of thin screens whose thickness are 
too small to support the intraslit resonance, the extraor- 
dinary transmission is caused by the excitation of sur- 
face plasmons or their matter-wave analog, surface mat- 
ter waves H, ^ IE 13' I'^ tl^is Letter, we show another 
mechanism that provides a great transmission enhance- 
ment not by coupling to the surface electron or matter 
waves but by interference of diffracted waves in the far 
field. 

The transmission enhancement without assistant of 
surface waves can be explained in terms of the follow- 
ing theoretical formulation. We first consider the trans- 
mission of light through a structure that is similar but 
simpler than an array of holes, namely an array of par- 



allel subwavelength slits. The structure consists of a 
thin metal screen with the slits separated by many wave- 
length. To exclude the plasmons from our model, the 
metal is considered to be a perfect conductor. Such a 
metal is described by the classic Drude model for which 
the plasmon frequency tends towards infinity. Owing to 
the great distance between the slits, the electromagnetic 
field at one slit is assumed to be independent from other 
slits. The transmission of the slit array is determined 
by calculating the light power in the far-field diffraction 
zone. The waves diffracted by each of the independent 
slits are found by using the Neerhoff and Mur approach, 
which uses a Green's function formalism for a rigorous 
numerical solution of Maxwell's equations for a single, 
isolated slit 0, 0, 0|. The calculations show a big, 
up to 5 times, resonant transmission enhancement near 
to the Fabry-Perot wavelengths determined by an array 
period. To clarify the numerical result, we then present 
an intuitively transparent analytical model, which quan- 
titatively explains the resonant enhancement in terms of 
the far-field interference of the propagating waves pro- 
duced by the independent shts. The model predicts the 
~5-times (resonant) and ^1000-times (nonresonant) en- 
hancements for both the light and matter waves passing 
through a perforated metallic or dielectric screen, inde- 
pendently on the apertures shape. Verification of the 
analytical formulae by comparison with data published 
in the literature supports these predictions. The Wood 
anomalies in transmission spectra of optical gratings, a 
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FIG. 1: The per-slit transmission Tm{X) of an array of inde- 
pendent slits of the period A versus the wavelength for dif- 
ferent number M of slits. There are three Fabry-Perot like 
resonances at the wavelenghts A„~A/n, n=l, 2 and 3. 



long standing problem in optics jlSj . follow naturally 
from interference properties of the model. 

Let us first investigate the light transmission by using 
the rigorous model. The model considers an array of M 
independent slits of width 2a and period A in a screen 
of thickness 6 <C A. The screen placed in vacuum is illu- 
minated by a normally incident TM-polarized wave with 
wavelength A = 2ttc/uj — 27r/fc. The magnetic field of the 
wave H{x,y,z,t) — U{x)exp{~i{kz + LL>t))ey is assumed 
to be time harmonic and constant in the y direction. The 
transmission of the slit array is determined by calculat- 
ing all the light power P(A) radiated into the far-field 
diffraction zone, xG[—oo, oo] at the distance z 3> A from 
the screen. The total per-slit transmission coefficient, 
which represents the per-slit enhancement in transmis- 
sion achieved by taking a single, isolated slit and placing 
it in an M-slit array, is then found by using an equation 
TmW = PW/MPi, where Pi is the power radiated by 
a single slit. Figure 1 shows the transmission coefficient 
TmW, h^ tfic spectral region 500-2000 nm, calculated for 
the array parameters: a = 100 nm, A = 1800 nm, and 
& = 5 X 10~^ Xmax- The transmitted power was computed 
by integrating the total energy fiux at the distance z — 1 
mm over the detector region of width Ax = 20 mm. The 
transmission spectra TmW is shown for different values 
of M. We notice that the spectra is periodically 

modulated, as a function of wavelength, below and above 
a level defined by the transmission ri(A) = 1 of one iso- 
lated slit. As M is increased from 2 to 10, the visibility of 
the modulation fringes increases approximately from 0.2 
to 0.7. The transmission Tm exhibits the Fabry-Perot 
like maxima around wavelengths A„ ^ A/n {n=l, 2, ...). 
The spectral peaks increase with increasing the number 
of slits and reach a saturation (T^"^ « 5) in amplitude 



by M = 300, at A « 1800 nm. The peak widths and the 
spectral shifts of the resonances from the Fabry-Perot 
wavelengths decrease with increasing the number M of 
slits. From the data of Fig. 1, one can understand that 
enhancement and suppression in the transmission spectra 
are the natural properties of an ensemble of independent 
subwavelength slits in a thin (6 <^ A) screen. The spec- 
tral peaks are characterized by asymmetric Fano-like pro- 
files. Such modulations in the transmission spectra are 
known as Wood's anomalies. The minima and maxima 
correspond to Rayleigh anomalies and Fano resonances, 
respectively 0| . A surprising prediction of the model is 
a weak Wood's anomaly, which is present in a classical 
Young type two-slit system (M = 2). 

The above-presented data is based on rigorous calcula- 
tion of the energy fiux by using the electromagnetic field 
evaluated numerically. The transmission enhancement is 
achieved by taking a single, isolated slit and placing it 
in an array. The interference of the waves diffracted by 
independent slits can be considered as a physical mech- 
anism responsible for the enhancement. To clarify the 
results of the computer code and gain physical insight 
into the enhancement mechanism, we have developed an 
analytical model, which yields simple formulae for the 
diffracted fields. For the fields diffracted by a narrow 
(2a ^ A, & > 0) slit into the region \z\ > 2a, it can 
be shown that the Neerhoff and Mur model simplifies to 
an analytical one. For the magnetic H = {0,Hy,0) and 
electric E — {E^, 0, E^) fields we found: 

Hy{x, z) = iaDF^{k[x'^ + z^J^/^), (1) 

E,{x,z) = -az[x^ + z^]-^/^DFl{k[x^ + z^]^/^), (2) 
and 

E,{x, z) = ax{x'' + 7?\-^l'^DFl{k{x' + z^Ji/Z)^ (3) 
where 

D = 4fc-i[[exp(ifcfe)(ayl - k)f - (aA + hfY^ (4) 

and 

A = Flika) + '^\h[ka)Fl(ka) + h(}ia)Fl(ka)\. (5) 

Here, F\, Fq, Fq and F\ are the Hankel and Struve 
functions, respectively. The fields are spatially nonuni- 
form, in contrast to a common opinion that a subwave- 
leng th aperture diffracts light in all directions uniformly 
[l6|. The fields produced by an array of M independent 
slits are given by E{x, z) — X^^^j^ Em{x + mA, z) and 
H{x, z) — X]m=i Hm{x + mA, z), where Em and Hm are 
the fields of an m-th beam generated by the respective 
slit. As an example. Fig. 2(a) compares the far-field dis- 
tributions calculated by the analytical formulae (1-5) to 
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FIG. 2: Electromagnetic fields in the far-field zone, (a) The 
fields Re{E^{x)) [A and D), RB{Hy{x)) {B and E), and 
Re{10Ez{x)) (C and F) calculated for M = 10 and A = 
1600 nm. The curves A, B, and C: rigorous model; curves 
D, E, and F: analytical model, (b) Re{Ex(x)) for M=\: 
analytical model, (c) Re{Ex(x)) for M=5: analytical model. 



that obtained by the rigorous model. We notice that the 
distributions are undistinguishable. Thus, the analytical 
model not only supports results of our rigorous model, 
but presents an intuitively transparent explanation of the 
enhancement in terms of the interference of the beams 
produced by the multi-beam source. The array-induced 
decrease of the central beam divergence (Figs. 2(b) and 
2(c)) is relevant to the beaming light ^-^id the non- 
diffractive light and matter beams [l8| . 

The analytical model accurately describes the fields 
E and H, but what is about the transmission coeffi- 
cient? The field power P, which determines the co- 
efficient Tm, is found by integrating the energy flux 
\S\ = \E X H* + E* X H\. Thus, the model accurately 
explains also the light transmission. We now consider the 
predictions in light of the key observations published in 
the literature for the two fundamental systems of wave 
optics, the single-slit and two-slit systems. The major 
features of a single-slit system are the intraslit resonances 
and the spectral shifts of the resonances from the Fabry- 
Perot wavelengths |^ . In agreement with the predictions 
0, the formula (4) shows that the transmission T = 
P/Po = (a/A:)[Re(D)]2 + [Im(D)]2 exhibits Fabry-Perot 
like maxima around wavelengths A„ = 2b/ n, where Pq is 
the power impinging on the slit opening. The enhance- 
ment and spectral shifts are explained by the wavelength 
dependent terms in the denominator of Eq. (4) . The en- 
hancement (T(Ai)Ri&/™ 0) 

is in contrast to the atten- 
uation predicted by the model . The Young type two- 
slit configuration is characterized by a sinusoidal mod- 
ulation of the transmission spectra T'2(A) 0,l23l- The 
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FIG. 3: The per-slit transmission coefficient T(X) versus 
wavelength for the Young type two-slit experiment |19| . Solid 
curve: experiment; dashed curve: analytical model. Parame- 
ters: a = 100 nm, A = 4900 nm, and b = 210 nm. 



modulation period is inversely proportional to the slit 
separation A. The visibility V of the fringes is of order 
0.2, independently of the slit separation. In our model, 
the transmission is given by T2^J[Fl{xi)[iFQ{xi)]* + 
Fl {X2)* iFq (x2)]dx , where xi = x and X2 = a; + A. The 
high-frequency modulations with the sideband-frequency 
/,(A) «/i(A) + /2(A,A)-1/A (Figs. 1 and 3) are pro- 
duced as in a classic heterodyne system by mixing two 
waves having different spatial frequencies, /i and /2. 
Although our model ignores the plasmons, its prediction 
for the visibility {V w 0.1) of the fringes and the resonant 
wavelengths \n = ^/n compare well with the plasmon- 
assisted Young's type experiment IQj (Fig. 3). In the 
case of 5 > A/2, the resonances at A„ = A/n can be 
accompanied by the intraslit resonances at A„ — 2b /n. 
Some difference between the calculated and measured 
values shows that the plasmonless and plasmon-assisted 
resonances can compete in certain situations. It can be 
mentioned that although our model considers a screen 
of perfect conductivity, polarization charges develop on 
the metal surface. The proponents of surface polari- 
tons, however, do not adhere strictly to traditional sur- 
face plasmons. We considered TM-polarized waves be- 
cause TE modes are cutoff by a thick slit. In the case 
of a thin screen, TE modes propagate into slits so that 
magneto-polaritons develop. Because of the symmetry of 
Maxwell's equations the scattering intensity is formally 
identical with E and H swapping roles. Again, no tradi- 
tional surface plasmons are present in our model. 

In order to gain physical insight into the mechanism 
of plasmonless enhancement in a multi-slit (M>2) sys- 
tem, we now consider the dependence of the transmis- 
sion Pa/ (A) on the slit separation A. We assume that 
the slits are independent also at A 0. According to 
the Van Citter-Zernike coherence theorem, a light source 
(even incoherent) of radius r = M{a + A) produces a 
transversally coherent wave at the distance z<7ri?r/A in 
the region of radius R. Thus, in the case of a subwave- 
length light source (A <C A), the collective emission of 
the ensemble of slits generates the coherent electric and 
magnetic fields, E = ^^^^^j^ Pmexp(i(y£>m)~Af Piexp(z(/?) 
and HthM HieyiY>{itf) . This means that the beams ar- 
rive at the detector with the same phases ipm{x) « v{x). 



4 




Wavelength ( nm ) 

FIG. 4: The per-slit transmission Tm{^) versus wavelength 
for the different values of A and M: (A) A = 100 nm, M = 2; 
(B) A = 500 nm, M = 2; (C) A = 3000 nm, M = 2; (D) A = 
100 nm, M = 5; (E) A = 500 nm, M = 5; (F) A = 3000 nm, 
M = 5. Parameters: a = 100 nm and 6 = 10 nm. There are 
two enhancement regimes at A <g; A and A>A. 



Consequently, the maximum power of the emitted hght 
scales with the square of the number of slits (beams), 
P - AP. Therefore, the transmission {Tm^P/M) grows 
linearly with the number of slits, Tm''^M. For a given 
M, the function Tm(A) monotonically varies with A. 
Such an enhancement regime (first regime^) is shown in 
Fig. 4. At appropriate conditions, the transmission can 
reach 1000-timcs enhancement (Af = Xz/nR^a + A)). 
The transmission enhancement in the second regime does 
not require close proximity of the slits. In the case of 
R > Xz/irr (A>A), the beams arrive at the detector with 
different phases Consequently, the power and trans- 
mission grow slowly with the number of slits (Figs. 1- 
4). According to our model, the transmission Tm ex- 
hibits the Fabry-Perot like maxima around wavelengths 
A„ = A/n. The constructive and destructive interfer- 
ence of the beams leads respectively to the enhancement 
and suppression of the transmission amplitudes as in a 
classical heterodyne system (Figs. 1 and 3). 

The analytical model gives not only intuitively trans- 
parent explanation of the plasmonless transmission en- 
hancement, but in contrast to the previous studies, pre- 
dicts the enhancement for the matter waves. Indeed, in 
our model, the enhancement is based on coherent ex- 
citation of an assemble of slits and interference of the 
diffracted waves (beams). The constructive interference 
is provided not by coupling between the slits, but by 
a geometrically well-defined phase relationship between 
wave amplitudes produced by a subwavelength source 
at different lateral locations in the far field. Thus, the 
enhancement mechanism depends neither on the nature 
(light or matter) of the fields Vm nor on material and 
shape of the apertures. For instance, in the first en- 
hancement regime, the fixed phase correlation leads to 



the M-times enhancement of the field amplitude, ijj = 
V'm exp(i(^TO)~-MV'i 6xp(i(p), and consequently to 
the M-time transmission enhancement. Due to Babinet's 
principle, the model predicts the enhancement also in the 
refiection spectra. Another surprising prediction of the 
model is a fact that the destructive interference of the 
fields at the detector can lead to the zero transmission. 
Indeed, the interference of the positive (-|-'0) and negative 
{—ip = V'exp(i(^+7r)) fields produces a field with the zero 
amplitude and energy. Two subwavelength glass fibers of 
different length can be used to generate such fields. The 
value T=0 is obtained by summing up the respective pos- 
itive and negative energies. Notice that the amplitudes 
of the fields ipm can rapidly decrease with increasing the 
distances a;, y and z. However, due to the enhancement 
and beaming mechanisms (Figs. 1-4), an array produces a 
propagating wave ip with low divergence. Such a behav- 
ior is in agreement with the Huygens-Presnel principle, 
which considers a propagating wave as a superposition of 
secondary spherical waves. 

The phenomena of enhanced and zero transmissions 
arc strange and paradoxical from a point of view of the 
energy conservation. Indeed, the energy conservation is 
a general law of nature. In the classical physics, the en- 
ergy Wab of a conservative system AB that contains 
the objects A and B is equal to a sum of the positive 
energies Wa and Wb of the isolated objects A and B, 
Wab = Wa + Wb- For an example, the kinetic energy 
Kab of the system of two macroscopic objects is given 
by Kab = p'a/^itia + p%/2mB- If the objects A and B 
arc two waves or they do have a wave nature of micro- 
scopic (quantum) objects, the energy conservation does 
not exhibit the law Wab = Wa + Wb- It can be eas- 
ily demonstrated that the energy of two harmonic light 
waves having the same amplitudes {Ea = Eb = Eq) and 
energies {Wa=Wb=Wo=EqEq^ dx) is given by Wab = 
Wa+Wb+IEqE^^ Jcos{ipa{x) - (pB{x))dx, where ipa{x) 
and <pb{x) are the wave phases at the coordinate x, and 
At = 1. The interference term, which is the additional 
energy from a point of view of the classical (non-wave) 
physics, is the positive or negative energy associated with 
the wave phases. In the case of two waves having the 
same phases {ipa{x) = (Pb{x)), the energy is given by 
Wab = Wo + Wo + 2Wo, where 2Wo is the "extra energf 
according to the classical physics. Correspondingly, the 
energy of the M waves is given by Wm = M'^Wq. An- 
other counter-intuitive feature of the two-wave system is 
the" annihilation of energ'i/^ {Wab = Wo+Wo+{—2Wo)) 
in the case of (Pa{x) — tpBix) = tt. Notice that the to- 
tal energy is always positive, while the "additional" en- 
ergy can be negative- The "negative energif' cannot be 
separated from the total positive energy of the system. 
Furthermore, the "negative energy" cannot be extracted 
from the system and transformed to the external sys- 
tem. Therefore, the " negative energy" can be considered 
as the internal virtual energy. Although, the construe- 
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tive and destructive interference are well known for the 
light intensity, up to now, the interference phenomenon 
for the energy has been considered as unphysical. We 
have showed the " non-classic energy conservation" based 
on the " extra energif^ , " negative energy" and " annihila- 
tion of energy" for multi-wave (multi-beam) systems in 
the context of extraordinary transmission of light and 
matter by subwavelength aperture arrays. 

It is worth noting that the presented model is similar 
in spirit to the dynamical Bloch-waves diffraction model 
[2l|. the Airy-like model based on the Rayleigh field ex- 
pansion [23, and especially to a diffracted evanescent 
wave model The difference arises from a fact that 

the models |llLl2llE3 | consider the case of M = 00, while 
predictions of our model strongly depend on the number 
M 7^ 00 of slits. In addition, our model deals with in- 
dependent slits, while the models of Refs. |^ 0, 
consider the slits electromagnetically coupled via the pe- 
riodic boundary conditions. There is an evident resem- 
blance also between our model and a Dicke superradi- 
ance model jl^l of collective emission of an ensemble of 
atoms. A quantum reformulation of our model can help 
us to understand why a quantum entangled state of pho- 
tons is preserved on passage through a hole array [2J|. 
The quantum model will be presented in our next paper. 
Notice, that the surface waves can couple the radiation 
phases of the different slits, so that they get synchronized, 
and a collective emission can release the stored energy 
as an enhanced radiation. This kind of enhancement is 
of different nature compared to our model because the 
model does not require coupling between the slits. 

In conclusion, we showed that the energy conserva- 
tion does not exhibit the "classic law" if the physical 
objects are waves or they do have a wave nature of micro- 
scopic (quantum) objects. The " extra energy" , " negative 
energy" and " annihilation of energy" are predicted for 
multi-wave (multi-beam) systems. The paradoxical phe- 
nomenon was demonstrated in context of the extraordi- 
nary transmission of light and matter through subwave- 
length aperture arrays assisted by surface waves [T.W. 
Ebbesen et al.. Nature (London) 391, 667 (1998)) and 
E. Moreno et al., Phys. Rev. Lett. 95, 170406 (2005)]. 
The physical mechanism of "non-classic energy conser- 
vation" provides the transmission enhancement (" extra 
energy" ) not by coupling to the surface waves but by in- 
terference between diffracted waves (beams) in the far 
field. Furthermore, the " annihilation of energy" by de- 
structive interference of fields at the detector leads to 
the zero transmission. Another surprising prediction of 
the model is a weak Wood's anomaly, which is present 
in a classical Young type two-slit system. The Wood 
anomalies in transmission spectra of optical gratings, a 
long standing problem in optics, follows naturally from 



interference properties of the model. We believe that the 
"non-classic energy conservation" demonstrated in the 
present paper gains physical inside into the nature of en- 
ergy conservation. The physical mechanism of the " extra 
energy" , " negative energy" and " annihilation of energy" 
affects the energy conservation in subwavelength systems. 
The mechanism can be used for numerous applications, 
for an example, in nanooptics and nanophotonics based 
on transmission and beaming of light or matter by sub- 
wavelengths apertures. 
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